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ON GENERALIZED SPACE DIFFERENTIATION OF 
THE SECOND ORDER. 

By B. O. Peirce. 

Presented December 9, 1903. Received December 26, 1903. 

If one has to investigate the strength of a field of force defined by a 
given scalar potential function, or to study the flow of electricity 
in a massive conductor under given conditions, or to apply Green's 
Theorem to given functions in the space bounded by a given closed 
surface, or, indeed, to treat any one of a large number of problems in 
Mathematical Physics or in Analysis, one often needs to find the numer- 
ical value at a point, of the derivative of a point function taken in a 
given direction. This has given rise to the familiar idea of simple space 
differentiation and of the normal derivative of one scalar function with 
respect to another ; indeed the properties of the first and of the higher 
space derivatives of a function of n variables taken with respect to any 
fixed direction in n dimensional space, have been treated very clearly 
and exhaustively by Czuber.* 

It is sometimes desirable to use also the conception of general space 
derivatives of the second order. This is the case, for instance, when 
one is determining the rate of change of the intensity of a conservative 
field of force at a point which is moving, either along a curved line of 
force or on a curved surface related to such lines in a prescribed man- 
ner. It is easy to define the general space derivative of any order of a 
given function. 

This paper discusses very briefly a few elementary facts with regard 
to generalized space differentiation of the second order, and treats first, 
for the sake of simplicity, differentiation of functions of two variables, 
in the plane of those variables. 

Plane Differentiation. 

Let there be in the xy plane two independent families of curves 
(u = c, v = h) such that in the domain, R, one and only one curve of 

* E. Czuber, Wienerberichte, p. 1417 (1892). 
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each family passes through every point and no curve of either family 
has anywhere a multiple point. At every point, P, in the domain, the 
two curves (one of the u family and one of the v family) which pass 
through the point indicate two directions, s 1} s 2 , and if the sense of each 
of these be determined by any convenient convention, they may be 
defined by pairs of direction cosines (J r , w^), (4, m^), where ly, »* 1} l 2 , m % 
are given scalar functions such that at every point 

If O is any scalar function of the codrdinates which within R has 
finite derivatives of the first and second orders with respect to these 
coordinates, the derivative of Q at P in the direction s x is the value 
at the point of the quantity 

h 'te + mi 9j (2) 

and this new scalar function of x and y may be conveniently indicated 
by the expression \D, H] P . If P' is a point on the u curve which passes 
through P, taken near P and in the sense of the direction * 1} [D s il\ p is 

l&pr — ilp 

the limit, as P' approaches P, of === 

If on the curve of the second (or v) family which passes through P, 
a point Q be taken near P and in the sense of the direction s 2 , the limit, 
as Q approaches P, of the quantity 

PQ (3) 

may be indicated by the expression \D, D, €t] P , and this is the second 
derivative of O at P taken with respect to the directions * t and s 2 in the 
order given. . 

Thus, if 

2 x 1 

O = 2 x 2 — y 2 , li= — ===== , m 1 = 



V^x^+l V*x*+1 

V^ + tf Vx* + y* ' V 4z 2 +l 

x (32 st? y + 8y 2 + 16y + 8x 2 + 2) 



A D.a 



Vx* + y 3 ■ (4a; 2 + 1)1 



PEIRCE. — SPACE DIFFERENTIATION OF THE SECOND ORDER. 379 
It is evident from the definition just given that 

9 2 il 9*0, 9' 2 Q 



+ 



( . 9h Sh\ 90 ( 9m, 9m, \ 90 



and that D„ D, O is quite different in general from D s D s O : the order 
of the two differentiations is material. 

If the u curves happen to be a family of parallel straight lines and the 
v curves another family of parallel straight lines, 

D H D h Q = kk-gz + (km, + km,) ■ g— ^ + Wl • Mj • _, (5) 

and the coefficients in this expression are constants. 

If the m curves and the v curves are identical and are a family of 
straight parallel lines, we have 

S 2 fi 5 2 3 2 fi 

A^ = V-^ + 2W.^ + ^.^, (6) 

the familiar form of the second derivative of O along the fixed direction 
«!, which often appears in work involving the transformation of Cartesian 
coordinates. Simple special cases of this formula are obtained by putting 
/equal to 1, 0, and m. 
Since I* + mf — 1 , 

9h 9m x 9m, 91% 

9x 9y 9x 9y ' 

90 
and if at any point s, and s 2 are such as to make the coefficient of yr- 

9o 

in (4) vanish, the coefficient of ^r— will vanish also. Such points as 

this lie, in general, on a definite curve, the equation of which is to be 
found by equating one of these coefficients to zero. If s t is a fixed di- 
rection so that I, and m x are constants, (4) takes the form (5), but the 
coefficients are not constants unless s 2 also is fixed. 

If the two variable directions s,, s 2 coincide, (4) becomes the second 
derivative of the function O taken with respect to the direction s, ; 
that is, 
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D li «0 = ( l -- s? +8l 1 . 1 . 5 _^. + w .._ 



+ 



(j 91, 9l\9a ( 9m, 9mA 9n .„ 

If the direction cosines of a plane curve at a point on it are I and m, 
the curvature of the curve at P has the same absolute value as have the 
expressions 

1 / 91 9l\ 1 / 3w , 3m\ 

»(/ fe + m *} 7(/"&- + ,B "*} < 8 > 

If, therefore, two directions, s u s s , are defined by two curves which, at 
a point, P, common to both, have a common tangent and equal curvatures, 
the second derivatives at P of a function Q taken with respect to the 
two directions are equal. 

If at any point the curvature of the curve of the w family which 
defines the direction s 1 is zero, the coefficients of 9Cl/9x and 90 /9y in 
the expression for D s 2 il at the point vanish. If the u curves are a 
family of straight Hues, the last two terms of (7) disappear, but the 
coefficients of the other terms are, in general, not constant. 

If there is no point in the region R at which both the quantities 3fl/ 9x, 
SQ \9y vanish together, and if the direction s is at every point of R that 
in which fi increases most rapidly, D S U = h, where h is the gradient of 
Q, that is, the tensor of the gradient vector. Now h is itself, in general, 
a scalar point function, which, when equated to a parameter, yields a fam- 
ily of curves the directions of which are usually quite different from those 
of the lines of the gradient-vector. The normal at any point P to the 
curve of this h family which passes through the point, has the direction 
cosines 



3x1 ' 9y 

where h! is the gradient of h. The angle between the direction, s, of the 
gradient vector of O and the normal to the h curve has at every point 
the value 

/o r> V 9h SO Sh 901 I ... 

and the second derivative of O with respect to the direction * is, therefore, 

«°-[s-5 + g-&]/*-*-<ft<>- < 10) 
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Let the normal derivative,* at auy point P, of a point function V, 
taken with respect to another point function W, be the limit, as P Q ap- 
proaches zero, of the ratio of Vq — V P to W Q — Wp, where Q is a point 
so chosen on the normal at P to the surface of constant W which passes 
through P, that W Q — W p is positive : if, then, ( V, W) denotes the angle 
between the directions in which V and W increase most rapidly, the 
normal derivatives of V with respect to W, and of W with respect to V, 
may be written 

[D w V]=h v - cob (V,W)/h„ [D 7 W~\ =h w - cos (V,W)/h v : (11) 

if h r = h m these derivatives are equal. 

With this notation (10) may be rewritten in the form 

D,*Q = h--[D Q h.] (12) 



If at any point D S 2 Q vanishes, it is easy to see from (10) that either 
the gradient (h 1 ) of h vanishes at the point, or else the h and O surfaces 
cut each other there orthogonally. This latter case is exemplified in the 
familiar instance of the electrostatic field due to two long parallel straight 
wires of the same diameter, charged to equal and opposite potentials : if 
the wires cut the xy plane normally at P t , P 2 , and if the line joining 
these intersections be taken for x axis with the point midway between 
them for origin, the potential function is of the form V= Alog^/r^, 
where r 2 = (x — af + y 2 , r 2 — (x + a) 2 + y 1 . The intensity of the 
field, in absolute value, is h = 2 a A/r l r 2 , and the second derivative of V 
taken along the line of force (that is, the rate at which the intensity of the 

field changes) is numerically equal to — -^ ^ ■ 

r x ' r% 

D 2 V taken along a line of force vanishes, therefore, at all points on 

they axis, and at all such points the curve of constant V (r l /r i = 6) 

cuts the curves of constant h (r x r 2 = k) orthogonally. At points on the 

y axis the direction of the lines of force is parallel to the x axis, and the 

second derivative of V with respect to the fixed direction x happens to 

0/ V 

vanish here also where I = 1, — = 0, m = 0, -^r— = 0. The quantity 

h! does not vanish at any finite point. 

* Peirce, The Newtonian Potential Function, p. 116. A Short Table of Inte- 
grals, p. 106. 
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The example just discussed is in contrast with the case where the O 
family are a set of parallel curves of any kind, and h in consequence (if 
not constant) is a function of O alone, so that the h curves and the O 
curves coincide, and if X>„ 2 Q vanishes anywhere, it must be where h! 
vanishes. A simple example of this is furnished by the field of attrac- 
tion within a very long cylinder of revolution, the density of which is a 
function of the distance from the axis alone. 

If the directions s t and s 2 are everywhere perpendicular to each other, 
we may without loss of generality write l 2 = — m x , m 2 = l l ; in which 
case the coefficients of SO/Sce, 9il/9y in (4) become 



v 2 9x 



9m 2 \ , /, 9l 2 9l 2 \ . 

+ w?2 ._j and _^._ + m2 ._j : (13) 

these vanish if the v curves form a family of straight lines, or the u 
curves a family of straight or curved parallels. The order of differentia- 
tion with respect to the orthogonal directions s ly s 2 is immaterial if both 
the u and the v curves are straight lines, that is, if the directions are 
fixed. 

If Sj is the direction in which Q increases most rapidly, and s 2 the 
direction of constant Q, 

n n n n /, P> n 9h 9il . 5h ~\ Ih 

i s*n f/soy fSayl 3a 9arm_9^f\ ) l h , U \ 
~ }9x^9yl\9x) \9y)_\ + 9x 9y \_9f 9x*]\l ' K } 

Now the direction cosines and the slope of the line of the gradient vector 
at any point are 

i 9a i 9a , 9ai9a 
h'~9i' h'dg' 9yJ9x- 

So that the curvature of the line is 



; _ r /5n/5oy"]i 

L + \Ssl9x) J 

( 3 2 o r/Soy /5oy-| 9ji.9nV9jn_9Hf\ j /., (1S) 
= 19x^1X9^) \9y) J + 9x 9y \_3tf 9**1 ) I y ' 
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and we may write in this case 

D, 2 D s Q = h/p. 

This expression gives the rate at which the maximum slope of the surface 
the coordinates of which are (x, y, O), changes as one goes along a line 
of level.* 

When s, and s 2 are perpendicular to each other, we have in general 

3 2 n s 2 n 3 2 o 

2 dx 1 dx • 9y 9y 

( 3m, 3m,\3n ( 91, 3l,\3n 

+ \ m ^te- h '^)te + { h '3y-- mi '9xj9y-' (16) 

and since l* + m, 2 = 1, 

, 3k . 3m, 3h , 3m, 

h • 7T- + m, ■ -=- = 0, I, ■ -s- + m, • ^— = 0. 
eta: 9x 9y 9y 

So that if we add together (7) and (16) we shall get 

Ai n + J o Si o__ + _ + _ ^ _ + _ _ _ 

2_ 3 2 Q g 2 Q 1 gm t 50 _ J^ % 3C1 

~ 3x 2 + 3y 2 I,' 3y ' 3x m,' Shc'dy' ( ' 

It is evident that the values of the space derivatives defined above are 
wholly independent of the particular system of rectangular coordinates 
which may be used. 

Space Differentiation. 

At every point of the space domain, R, let two independent directions 
(s„ s 2 ) be defined by the direction cosines (I,, m„ n,), (l 2 , m 2 , h 2 ), where 
l„ m„ n,, 4, in* «2 are any six single-valued point functions which 
satisfy the identities 

h* + m, 2 + V = 1, Z 2 2 + m 2 2 + w 2 2 = 1, (17) 

and have finite derivatives of the first order with respect to the coordi- 
nates x, y, z. If, then, O is any single-valued function of the coordinates 
which within R has finite derivatives of the first and second orders with 

* Boussinesq, Cours d'Analyse Infinitesimale, T. 1, f. 2, p. 236. 
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respect to these coordinates, the derivative of Q at the point P, in the 
direction s u is the value at P of the quantity 

n n _, 9a 9a 9a - 

Through the point P passes a curve of the Lmily defined by the 
equations 

dx = dy = dz^ . 

l 2 m 2 ra 2 ' 

and this curve indicates the direction s 2 . If on this curve a point <2 be 
taken near P and in the sense of the direction s 2 , the limit, as Q ap- 
proaches P, of the quantity 

[AOWA.O], . (20 ) 

may be represented by [Z)„ 2?, fi] p and this is the second directional 
derivative at P of a taken with respect to the directions s t and s 2 in the 
order given. It is evident that 

\ 3a; cty 3z / c^ 

/, 3m, 5mx , 9m 1 \ 9a 

+ V 4 '^ + " 2 5F + % 5^J^' (21) 

and that this is not equal to D s D s O. 

If the directions s v s 2 are fixed, the six direction cosines are constants, 
the last three terms of (21) disappear, and the coefficients of the other 
six terms are constant. If the fixed directions * ls s 2 coincide, (21) re- 
duces to the familiar form 
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nso ,, 5 2 o 2 3 2 o , t 3 2 o 

+ 2/ lMl -^L-+2m 1% -^ r+ 2Z 1 n 1 -^ r , (22) 
ox • dy dy dz dx • dz ' 

whereas, if s 1 is not fixed, 

nao 72 3 2 n , 2 3 2 o , „ S 2 o , ol 3 2 o . - S 2 o 

+ 2ZlWl ftP* + I* " * + " I dy- + ni Fz)9x- 
, /, Smj Swi , SmA SO 

+ ( v — + m i ' ^vr + % - -~— I --- 



Sx Sy dz J 9y 



( I •— 1 

V 1 Sx 



S»i SwASO . 

+ Wl s? + ni s^)s^ (23 > 



All the coefficients in (22) are constants ; all those of (23) are in 
general variable. If s t is defined by any infinite system of straight lines 
of which just one passes through every point of space, and if the direc- 
tion «j at all points of any one of the lines is that of the line itself, the 
coefficients of SO / Sx, SO / dy, 30 j dz in (23) vanish. In particular, 
if the direction s 1 is that of the radius vector from a fixed point (a, b, e), 
(23) takes the form of (22) though the remaining coefficients are not 
constants. In any case if the coefficients of two of the three quantities 
30 / Sx, 30 / 3y, 30 / dz vanish, the third must vanish also. 

If the gradient, h, of O does not vanish at any point of R and if s is 
the direction in which O increases most rapidly, 

D,n = h, 

(24) 
£>20- rS^ . SO SA _ SO 3A _ SO"] / , 

[_Sx Sx 3y 9y dz dz J / 

If h! is the gradient of the scalar point function which gives the value 
of h, and if (O, h) represents the angle between the directions in which 
the point functions O and h increase most rapidly, 

/o « Y 3h 5n , 5h 50 , dh SOI / . „ 

and Z> s 2 O = h' ■ cos (O, h) , or h [D a h] (26) 

vol. xxxix. — 25 
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where \_Dah'] represents the normal derivative of h with respect 
to O. 

If the equation CI = k happens to represent a set of parallel surfaces, 
h, if not constant, is a function of Q alone, so that the h and surfaces 
are coincident : in this case cos (O, h) = 1 and D„ 2 O can vanish only 
where h! vanishes. In general, D? £2 vanishes when the h and Q sur- 
faces cut each other at right angles- 

If «!, « 2 ) *a ar e any three mutually perpendicular directions, 
h 2 + V + h 2 — m i 2 + m 2 + w s 2 = «i 2 + «2 2 + n s 2 = 1, 
h m x + l 2 »j 2 + l 3 m z = m x n x + m i n i + m 3 n s — l x n l + l 2 n 2 + l 3 n s = , 

9 2 n 3 2 o 3 2 o 

and Z yo + ZVO + D % «0=! ? + ^ + ! ? 

, SO r 3/, , 34 , 3Z 3 , 9l x ^ $k ^ 91,1 

dQ, r 9m x 9m 2 9m s 3m 1 7 9m 2 , , 9m s l 

, 30 r, 3w x , , 3» 2 , , 3w 3 , 3% , 3% 3n 3 H 

+ ^'9^ +k -^ +k '^ +m ^w +m2 'w +nh "^r 

(27) 
The Jefferson Physical Laboratory. 



